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The Arak process is a solvable stochastic process which generates coloured pat- 
terns in the plane. Patterns are made up of a variable number of random non- 
intersecting polygons. We show that the distribution of Arak process states is 
the Gibbs distribution of its states in thermodynamic equilibrium in the grand 
canonical ensemble. The sequence of Gibbs distributions form a new model pa- 
rameterised by temperature. We prove that there is a phase transition in this 
model, for some non-zero temperature. We illustrate this conclusion with sim- 
ulation results. We measure the critical exponents of this off-lattice model and 
find they are consistent with those of the Ising model in two dimensions. 



Figure Captions 



Figure |l| (A) A state \ of the Arak process (B) The discontinuity set 7 of (A). 

Figure |2| (A) A set of lines I intersecting T> (B) an admissible graph drawn on 
the set I (C) one of the two colourings of T> with discontinuity set given 
by the graph in (B). 

Figure |3| Updates in the Markov Chain Monte Carlo. Dashed and solid edges 
are exchanged by the moves, which are reversible. (A) Interior vertex birth 
and death (B) move a vertex, and (C) recolour a region by swapping a pair 
of edges. In an extra move, not shown, a small triangle may be created or 
deleted. Further move types are used to update boundary structures. 

Figure ^ Binder parameter Ud (see text), regressed with cubic polynomials. 
Curves correspond to distinct box-side lengths d. The maximum likelihood 
fit, constrained to intersect at a point, is shown. Error bars in this and all 
other graphs are la. 

Figure |js] The Binder parameter data of Figure [| rescaled with Ising critical 
exponents. The regression is a cubic polynomial. X43-4 — 38.5 for the fit 
is acceptable. 

Figure || The magnetisation m,d(T), regressed with cubic polynomials. 

Figure |7| The magnetisation data of Figure || rescaled with Ising critical ex- 
ponents. The regression is a quartic polynomial. The value of the \ 2 
statistic shows that the fit is a poor one. 

Figure ||| A selection of states equilibrated in a box of side d — 12 at tempera- 
tures below and above the estimated critical temperature T c ~ 0.6665(5). 



1 Introduction 



The Widom-Rowlinson model, with two species of discs and hard-core interac- 
tions between discs of unlike species, is sometimes referred to as the "continuum 
Ising model" . However there is another continuum model which might share the 
title. In 1982 Arak |Q presented a stochastic process in the plane with realisa- 
tions of the kind shown in Figure |l]A. States are composed of a variable number 
of coloured non-intersecting random polygons. Remarkably, the normalising 
constant is available as an explicit function of the area and boundary length of 
the region in which the process is realised. We present rigorous results and sim- 
ulation based measurements related to critical phenomena in a two dimensional 
"continuum Ising model" derived from the Arak process. 

There are few rigorous results for continuum models of interacting extended 
two dimensional objects. Moreover, relatively few Monte Carlo simulation stud- 
ies have been made, perhaps on account of the complexity of the simulation 
algorithms required. The Widom-Rowlinson model has a phase transition Q. 
Its critical exponents have been measured and put it in the Ising universality 
class Critical phenomena are known to occur in a range of related mod- 
els with q > 2 species and certain soft-core interactions j|, Where critical 
exponents have been measured |6| the universality class seems to be the class 
of the corresponding q~ species Potts model. For single-species models rigorous 
existence results for phase transitions have been given only in certain restricted 
models having area interactions ]?], || . 

In the model we consider the interface between black and white regions 
summarises the state in the same way that Peierls' contours parameterise an 
Ising system. The energy associated with a state is proportional to the length 
of the interface. In contrast to the Ising model, the vertices of the polygon 
forming the interface take positions in the plane continuum. At a temperature 
T = 1, the model we consider corresponds to the Arak process. For this value 
of the temperature the partition function equals the normalising constant of the 
corresponding Arak process. At smaller values of the temperature we are no 
longer dealing with an Arak process. We no longer have a closed form for the 
partition function. However the model remains well defined, and two phases 
coexist at temperatures bounded away from zero. 

Besides this result, which we prove, we estimate the critical exponents of the 
temperature-modified Arak model, using Markov chain simulation to generate 
realisations of the process. Values (obtained by "data-collapsing") are consistent 
with the corresponding critical exponents of the Ising model. This is in accord 
with what we expect from the hypothesis of universality, since the ground state 
of the temperature-modified Arak model is two-fold degenerate, and states are 
two dimensional. 

Although there is no high temperature limit for polygonal models (a class of 
models including the Arak process) consistent polygonal models might play this 
role (this point is made in ||). We give no rigorously determined upper bound 
on the critical temperature, although it is clear, from our simulations, that the 
consistent Arak process has a single phase. 
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2 The Arak process 



We now define the Arak process, following |l0| |. A state is a colouring map 
X ■ T> — » J from each point in an open convex set V C Sft 2 , onto a set J 
of possible colours. See Figure [l]A. We write dT> for the set of points in the 
boundary of V. We consider the simplest case, J — {black, white}, of two 
colours. 

Let Xx> be the class of all finite subsets x of T> U dV. Let (n > 1) 

(n) 

be the set of point-sets x composed of n points, so that X-p = U^L Xp , with 
= {(/)} the subset x = <f> containing no points. Let dxi be the element of 
area in T> and length on dT>. A measure dv(x) is defined on Xd by 

dv(x) — dx\dx2 ■ ■ ■ dx n . 

This is the measure of an independent pair of Poisson point processes of unit 
intensity, on the boundary and interior. 

Let Tx>(x) be the set of all "polygon graphs" 7 which can be drawn on the 
point-set x, ie the set of all graphs which can be drawn in T> with edges non- 
intersecting straight lines, with the points in x as vertices. All interior vertices 
must have degree 2 (they are V- vertices), and all boundary vertices degree 1 
(/-vertices). 7 is composed of a number of separate polygons which may be 
chopped off by the boundary. See Figure [l]B. 

The space of all allowed polygon graphs is the union over vertex sets x of 
the polygon graphs of x: 

r v = |J r v (x). 

xEX-d 

We define a measure on Tt> by 

dA( 7 ) = k( 7 )<Mz(7)), (!) 




for a pattern 7 with vertices at ^(7) = (x\, x% . . . x n ). is the smaller angle at 
vertex i in 2?, and the smaller angle made with the boundary tangent at x\ for 
vertices on dT>. The product over <i,j> runs over vertex pairs i,j connected 
by an edge in 7, with = \xi — Xj\ the length of the edge between vertices 
i and j. A counting measure is taken on the graphs of a fixed point set. The 
significance of k is sketched at the end of this section. 
Arak's probability measure on Tx> is 

Pvid-i) = ±- exp(-2L( 7 )) dA( 7 ), (3) 

with L(j) the summed length of all edges in 7, and 2-p a normalising constant. 
Remarkably, Zt> has a simple closed form [w| (ie the model is solvable), 

Z v = exp(L(dV) + irA(V)), 
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where L(dT>) and A(T>) are respectively the perimeter length and area of T>. 
Certain expectation values have been calculated (see [[l0| 0)- Some examples 
are given m Table | 

A colouring map \ : T> — > J is a function assigning a colour, black or white, 
to each point in T>. See Figure |l|A. Let a colouring map x be given and let 
B x be the set of points x G T> with a black point, ie some y 6 T> such that 
x(y) = black, in every e-neighbourhood. Let W x be similarly defined for white 
points. Let ~f(x) = B x D W x denote the discontinuity set of this colouring. For 
each polygon graph 7 we consider two colouring maps x '■ D — > J each having 
discontinuity set ■y(x) — 7- The two distinct colourings of a given polygon graph 
are assigned equal probability, so the probability measure for colour maps is just 
P v {d 7 ( X )}/2. 

The probability measure (||) has a number of beautiful properties, besides 
solvability. Striking are consistency and the Markov property. Consider an 
open region S of T> with e-neighbourhood (S% C T>; the probability measure 
for events in S, given full information about x on {dS) e , is independent of any 
further information about the state in V \ S. That is the Markov property. 
Next, let S be an open convex region S C T> and let 75 € L5 denote the 
restriction of a state 7 e Tx> to S. The probability measure for events simulated 
in V from Pv{dj} but observed in the subset 5" is equal to Ps{e?7}, in other 
words Psidj} = Pv{d~fs}- That is consistency. The Arak process shares 
these properties with a much larger family of probability measures called the 
consistent polygonal models. See for the general picture. 

We will now explain in brief how ^(7) arises, following JlO| closely. Consider 
a number of straight lines drawn in the plane. Let li = (pi, <pi) where pi is the 
perpendicular distance from the line to an origin and 4>i is the angle the line 
makes to the x-axis. The parameter space of U is L = [0, 00) x [0,7r). Let 
be that subset of L consisting of all lines intersecting V. Let dl = dp dip be 
Lebesgue measure of Lo- Let L 7 ^ be the set of all line sets £ = {h,h • ■ - In} 
made up of n lines, each in Lb- In this parameterisation Ct> = U n L D is the set 
of all sets of lines in the plane intersecting T>, and 

dv(t) = dl\dl2 ■ ■ ■ dl n 

is the element of measure of a line process in T>, corresponding to a Poisson point 
process of unit intensity in Ljj . Referring to Figure |5[ we define an admissible 
graph on a line set £ to be a graph with edges coinciding with lines in £, such 
that each line in £ contributes a single closed segment of non-zero length to the 
graph. All interior vertices are V vertices, all boundary vertices are / vertices. 
The set of all admissible graphs which can be drawn on some line set in Cx> is 
identical to T-p. Let 7 be some legal graph drawn on the line set £. Define a 
measure d\("f) — dv{£) in T-p using the line process as our base measure, and 
taking counting measure over the legal graphs of a line set. We now have two 
parameterisations of the graph: from its line set £, or from its vertex set x. The 
authors of have shown that d\(j) = dX^), ie, ^(7) arises as the Jacobian 
of the transformation between x and £. 
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3 Properties of a temperature modified Arak 



process 

We choose to modify the measure (|3|), and consequently loose solvability. Con- 
sider a system of non-overlapping polygonal chains of fluctuating number, length 
and vertex composition, confined to a planar region T>. The chains may be at- 
tached in some places to the boundary of V. The state is described by a graph 
7 G Yx>- Micro-states are associated with elements of volume d\("f) in Tp, so 
that in the Gibbs ensemble edge segments are isotropic in orientation (a rather 
unnatural choice). However, the Gibbs distribution Qv{dj} of this system is 
just the Arak distribution above, modified by the addition of a temperature 
parameter, as we now show. 

The Gibbs distribution Qvidj} has a density, 3(7) say, with respect to 
dA(7), the line measure. The Boltzmann entropy of the system is 

H[g}= [ g(i)Hg(i))d\(i) 

JT-D 

In the grand canonical ensemble, the energy and dimension of the system state 
fluctuate about fixed average values. We suppose that the state energy £'(7) is 
given by the total length of the chains, E(pf) — cL( 7 ), with c a positive constant. 
The dimension of the vertex position vector x is dim(x) = 2rii + with n$ (rib) 
the number of interior (boundary) vertices in 7. Maximising the entropy subject 
to constraints on the mean energy and mean dimension of the state, we obtain 
the distribution of systems of chains, 

Qv{<h} = ^cxp(- C L( 7 )/T) q- n °d\( 7 ), 

£>T>T 

where T and q are Lagrange multipliers, and n e is the number of edges in 7 
(n e = Hi + nt,/ 2). Under the change of scale Xi — > qxi, the measure transforms 
as dA( 7 ) — > q n 'd\(~f). We therefore set q = 1 without loss of generality. Setting 
c = 2 we obtain a "temperature-modified" Arak process 

P T ,v{dri = J— exp(-2L( 7 )/T) dA( 7 ). (4) 

-OT>T 

The function 2L( 7 )/T is a potential, (ie Z,x> T is finite), at least when < T < 1, 
and, by Theorem 8.1 of [pL2| , the temperature-modified measure keeps the spatial 
Markov property of the Arak measure. 

Let /i£(T) be the mean proportion of V coloured black (and /i^(T) white), 

fi%(T) = E T , v {A(B x )/A(D)}. 

The magnetisation of a state 

m( X ) = \A(B x )-A(W x )\/A(D) 

measures the colour asymmetry in that state. In our simulations (reported 
below) we see a qualitatively Ising-like temperature dependence in the mean 
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magnetisation. We prove, in an Appendix, that there is long range order (ie 
phase coexistence) in magnetisation, at all sufficiently small temperatures. We 
have translated Griffiths' version |D| of Peierls' proof of phase coexistence in 
the Ising lattice model to this continuum case. 

Let fjLp W (T) be the expected proportion of V coloured black given that the 
boundary is white, that is 

nT(T) = Et.v {A{B x )/A{D)\dD n B x = 0} . 

Theorem For the temperature modified Arak process in an open convex region 
T> C 5ft 2 there exists a temperature T C oid; < T co \^ < 1 and a constant a, a > 0, 
such that 

B\W frp \ s i 

\ 1 cold ) S 2 — a 

independent of the area A(T>) of the region. 

Surgailis || has shown that, for an open convex set S C T>, the thermo- 
dynamic limit V y 5ft 2 of Pr,r>{d7s} exists, for a class of measures including 
Pr,T>{dr/}, for all temperatures below some small fixed positive value. With the 
theorem above, 

B\Wlrp\ 1- B\W /rpx 

(i (T) = lim [i D ' (T) 

exists and satisfies /i B ' M/ (T co id) < 1/2 — a for some a > 0. Hence, there is phase 
coexistence at all temperatures T < T co \a. 

In fact it follows from the result stated in the Appendix that 

"""m^ (?+?+§)■ a» 

where z = (1/(tvT) — 1). Sketching the function of T on the right hand side of 
Equation @, we see that T co id > 0.18, though this bound is not at all sharp. 
Simulation (see below) shows that the model has a phase transition with critical 
temperature very close to T = 2/3. 

The proof of the theorem is in two parts. We are after an upper bound 
on the expected area coloured black. The area of black in a state with white 
boundaries is not more than the summed area of the polygons it contains, and 
is maximised when they are not nested. This observation leads to a simplified 
bound on the expected area coloured black, Equation (|l^). This first result is 
obtained by an obvious translation of the Griffiths calculation into the terms of 
a continuum process. In that case the next step, bounding the number ways a 
polygon can be drawn on a lattice of fixed size, using a fixed number of links, 
is straightforward. In the continuum, the analogous problem is to bound the 
volume of the parameter space of a polygon of fixed length, where volume is 
measured using A, the line-based measure. The main difficulty lies in the fact 
that there are unbounded, but integrable, functions in the measure which arise, 
for example, when an edge length goes to zero; these would be absent if there 
were no polygon closure constraint; as a consequence the closure constraint may 
not be relaxed as simply as it is in the Ising case. 
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4 Simulation Results 



The probability measure frz>{^7} may be sampled using the Metropolis-Hastings 
algorithm, and Markov chain Monte Carlo. In our simulations we take V to be a 
square box of side length d. Note that the number of vertices is not fixed. Since 
the dimension of a state depends on the number of vertices in it, the Markov 
chain must make jumps, corresponding to vertex addition and deletion, between 
states of unequal dimension. Simulation algorithms of this kind are widely used 
in physical chemistry Q, |1| and statistics [ijj]. Although there exist ver- 
tex birth and death moves sufficient for ergodicity, we allow a number of other 
moves in order to reduce the correlation time of the chain. See Figure |^. At 
each update we generate a candidate state 7', by selecting one of the moves, 
and applying it to a randomly selected part of the graph. The candidate state 
becomes the current state (ie it is accepted) with a probability given by the 
Metropolis-Hastings prescription. Otherwise it is rejected and the current state 
is not changed. In this way a reversible Markov chain is simulated. The chain 
is ergodic, with equilibrium measure Pr,v{dl}- Full details of our algorithm, 
including explicit detailed balance calculations for all the Markov chain updates, 
are given in [jn| . 

The sampling algorithm is quite complex, but because the model is solvable 
at T = 1, it is possible to debug the code, by comparing a range of estimated 
expectations with predicted values. In Table [l] we present a selection of system 
statistics at T = 1. Quantities in brackets are one standard deviation in the 
place of the last quoted digit. These data show how the simulation code was 
tested in comparisons using analytically derived expectation values. Let / equal 
the average of some statistic f(x) over an output sequence of length N, let p/(t) 
equal the normalised autocorrelation (or ACF) of / at lag t and, for M > 0, let 
Tf = 1 + 2 Y^tLi P(t) estimate the normalised autocorrelation time of /(x) in the 
output, so that the variance of / is estimated by r/var(/)/JV. We used Geyer's 
initial monotone indicator [Q to determine M, the lag at which the ACF is 
truncated. The asymptotic variance a 2 p of the ACF as t — > 00 was estimated and 
used as a consistency check on each measurement: the estimated ACF should 
fall off to zero smoothly, and at large lag should stay within 2a p bounds of zero. 
As usual we cannot show the Markov simulation process has converged, but it 
is at least stationary. 

Run parameters for the measurements at T < 1 are summarised in Ta- 
ble ||. Autocorrelations reported are for T = 0.66, near the critical temperature. 
We estimate the integrated autocorrelation time T m of the state magnetisation, 
along with its standard error jl9| and present these alongside the total run 
length. Referring to Table |^, the autocorrelation time is fitted within standard 
error by r m oc d 4 6 . Our Metropolis Hastings algorithm is a local update algo- 
rithm and this places practical limits on the size of the largest system we can 
explore. 

We now report our measurements of the mean magnetisation, fhd(T) = 
ET,d{m(x)}i an d the Binder parameter 

U d {T) = l- E ^ TO W 4 > 



3E T , d {m(x) 2 } 2 ' 
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Hi) 






Et= M {/(7)} 


ird 


Ad + 47rd 2 


4nd 2 




d 


L 


h e 


hi 


0.5 


1.571(6) 


5.13(2) 


3.13(2) 


1 


3.14(1) 


16.46(7) 


12.47(6) 


2 


6.27(2) 


58.1(3) 


50.1(2) 


4 


12.55(2) 


216.7(4) 


200.6(4) 


8 


25.14(2) 


836.2(5) 


804.2(5) 


xl 


1.1 


4.0 


5.1 



Table 1 : Listed are a selection of estimates made from output at T = 1 . Here d 
is the box side, and for a state 7, L(j) is the total edge length, n e (j) equals the 
number of edges, and 71,(7) equals the number of interior vertices. Quantities in 
brackets are one standard deviation and in the place of the last quoted digits. 

Under the scaling hypothesis, the various curves Ud{T) indexed by d all intersect 
at a single T- value, the critical temperature [po| , T — T c say. A Bayesian 
estimate T c may be given for the intersection point. Let U denote the ordered 
set of independent [/-measurements we made (43 in all), let Tjj denote the 
ordered set of T values at which measurements were made, and let denote 
the ordered set of estimated standard errors for the measurements in U. These 
data are represented by the error bars in Figure ^. Each measurement is an 
independent measurement. For each d = 6, 8, 12, 16, we model the unknown 
true curve Ud(T) using a cubic 

2 

U* d (T) = U* + (T-T*)J24 d)xP - 

The parameterisation constrains the regression in such a way that the four curves 
intersect at a point (T*, U*). We simulate the joint posterior distribution of the 
random variables 

a< 6 > , a ( 8 > , a( 12 ) , a( 16 > , T* , U* I U, T v , , 

conditioning the slope to be negative in the region containing the data, and 
conditioning the lines to intersect at a point, but otherwise taking an improper 
prior equal to a constant for all vectors of parameter values. Again MCMC 
simulation was used. The marginal posterior distribution of T* is very nearly 
Gaussian. Our estimate of the critical temperature is then 

f c = 0.6665(5). 

The quoted standard error is the standard deviation of T* in its marginal pos- 
terior distribution. 
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d 


# Updates 


7~m 




X 1U 




1 


16 


0.00181(7) 


2 


40 


0.018(1) 


4 


6400 


0.41(1) 


6 


16000 


2.4(3) 


8 


128000 


9.2(2) 


12 


300000 


59(3) 


16 


300000 


240(40) 



Table 2: Listed are run parameters for simulations at T = 0.66, a temperature 
close to the measured critical temperature. An update is a single pass through 
the Metropolis-Hastings propose/ accept simulation sequence. Measurements 
made at the same d value, but different temperatures, are based on the same 
number of updates. 

The Bayesian inference scheme used to estimate T c above is attractive for 
several reasons. Above all it quantifies the uncertainty in our estimate of T c , 
taking full account of the complex constraints applying in the regression (though 
taking no account of possible errors due to violations of scaling). The sensitiv- 
ity of the outcome to the orders of the regressing polynomials was explored. 
The chosen orders were the smallest that gave an acceptable likelihood. The 
posterior mode, which is the maximum likelihood estimate for T c , on account 
of our flat prior, occurs at T* = 0.6663. Metric factors weight the mass of 
probability in the posterior distribution only slightly away from the maximum 
of the likelihood. 

Because the energy has a discrete two-fold symmetry, and states are two 
dimensional, we expect the model to lie in the universality class of the Ising 
model. Finite size scaling under the scaling hypothesis leads to a system size 
dependence of the form |2(J 

fh d (r) = d-P^gid 1 /^) 
U d (r) = f(d^r) 

with / and g unknown functions, r the reduced temperature (T/T c — 1), and 
(3 and v critical exponents. If we plot Udij) or d^^rhdir) against d 1 / l/ r, we 
expect to see no significant dependence on system size d for r near zero. Using 
the Ising critical exponents v = 1 and (3 — 0.125 and our estimate T c for the 
critical temperature, we show, in Figures |^ and ^ the maximum likelihood fit 
to the transformed data. The transformed I/^-data lies on a smooth curve. The 
transformed m^-data does not give a satisfactory \ 2 (all of the misfit comes from 
points at T > T c ), but this is to be expected: we are seeing scaling violations (a 
satisfactory fit to a quartic can be obtained (xig-5 = 30) by dropping points at 
large T from the d — 6 and d — 8 data). If this is so, then the critical exponents 
of the Ising model the temperature dependent Arak process are equal at the 
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precision of our simulation analysis. 

Sample realisations from the model, taken at temperatures around the crit- 
ical temperature are shown in Figure |[ 
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Appendix: long range order 

We now give the proof. Condition on a white boundary. There can be no 
boundary vertices. Let Ty/v be the subspace of T-p of polygon graphs with no 
boundary vertices. Let 0x> be the subspace of Twd of graphs made up of just 
one polygon. Each point in 0£> corresponds to a single polygon, lying wholly 
in T>. We begin by proving the inequality Equation ( |To| ) below. 

Among states built from a given set of polygons, with no edge connected 
to the boundary, the black area is largest when the polygons are arranged so 
that none are nested. It follows that the area of black in a state \ with a white 
boundary is less than or equal to the sum of the areas of all the polygons in 
that state. The area of a polygon 9 of perimeter length L(9) is smaller than the 
area of a circle with the same perimeter, so A(9) < L(9) 2 /Air and 

MB X ) < E ^ W 

ec 7 (x) 

We want to take expectations of either side of Equation @ so we clear 7 from 
the domain of the sum, using 

£ f(9) = I f(9) 6(9 c 7 (X» d»{x(B)). 
ec 7 (x) 



Hi 



5(9 C 7(x)) puts a delta function at each point in corresponding to a polygon 
in 7. Each of these is a product of delta functions in T> for the vertices of 9 to 
coincide with those of 7, with an indicator function for the edge connections to 
coincide. x(9) is the set of vertex coordinate variables of the polygon 9. 

Now take the expectation of A(B X ) / A(V) over patterns x m ^wv- We have 

P>T < J q E{ 8(9 c 7 (X)) \dVC\B x =®} d v (x(9)). 

The expectation of the delta function is by definition 

L 5(9 C 7) x e ~ 2L ^/ T dXh) 
E{5(6c 7 )\dVnB x = n= J ^ J 2m/T . (7) 

Simplify the denominator by restricting the integral to those graphs to which 
the polygon 9 could be added without intersecting an edge of a polygon already 
in place. That is, if 

T 9 WV = {7 € T wv : 7 D 9} 
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is the set of polygon graphs containing the polygon 9, then 
is the sub-domain of interest. We have 

e -2L( 7 )/T dA(7) > I e -2L( 7 )/T dA(7) (8) 



We now turn to the numerator of Equation (Q). Carrying out the integration 
over vertices in 9 using the ^-function, 



3(9 C 7) x e~ 2L ^ )/T dX(-f) = / e- 2L ^/ T K{9) dX(-y\ 9) 

W<D JV wt> 

= n{9)e- 2L ^/ T ( e - 2i ^/ T dA( 7 ), (9) 



since 7 does not contain 9 in the second line. Substituting with (||) and (^|) in 
(0), and cancelling, 

E{ 5(9 C 7) I dV n B x = } < «(0) e - 2i ( )/ T , 

and consequently, 

B ^ <—L—f L(9) 2 e- 2L ^/ T dX(9). 



^ ~AnA(V)J e ; 
In close analogy with Griffiths' proof, we obtain 



B\W < 1 f u2 -2b/T 

^ ~ 4nA(V)J € 



S(b - L{6)) d\{9) 



db (10) 



The integral over b is an integral over polygon perimeter lengths. The problem 
is now to bound the integral over 0-p without introducing more than one factor 
of A(T>), or too rapidly growing a function of b. This is done by the following 

(n) 

Lemma. Let 0^-, be the subset of 0p of polygons with n vertices. 
Lemma Let 



so that 



m The 



and consequently 



J n = 5(b- L(9)) dX(6) (11) 
/ 5(b - L(6)) d\(fi) = J2 J„ 

J &v n=3 

(n—\n—ljjn—3 

J n <A{V)n 2 {n-l) ^ n J 2)] , 
S(b - L{9)) dX(9) < (2ir) 2 A(V)(4 + 2nb) 2 e 27rb . (12) 
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Proof of the Lemma: Start with J„ defined in Equation (fTl]). Use a 
standard labelling with x\ the variable corresponding to the vertex in 9 with 
the smallest x-coordinate, (smallest y-coordinate in case of ties) and vertex 
number increasing clockwise around 9. In the first step we break the polygon 
at x\ to make a chain. Consider the set O^ of distinct non-intersecting chains 
9 of n edges linking n+ 1 vertices, labeled with variables X\ to x n +i- All the 
vertices in a chain lie entirely to the right of the first vertex (or directly above). 
Polygons are chains, Op C Op , since the first and last vertices in a chain 
may coincide. Transform variables from {xi}™ =1 to {x\, {ej}^}, where is 
a Cartesian vector with origin Xi corresponding to the edge from the i'th to 
the (i + l)'th vertex. When we switch to integrating over chains, we constrain 
fil +£2 + • • • + £n to be zero, so that the polygon closes. Equation ( |ll| ) becomes 

J n < f S(b L{6)) 6^(Z k e k ) rf - ld - 2 --- d ~" dx u 
JqM eie 2 ...e„ 

with ei = L{e i ) and using sin(^j) < 1. 

The integrand is unbounded. We partition the space into regions, and impose 
the constraints b = L(9) and e Y + e 2 + . . . + e n = by integration over different 
variables in each region. For any particular region, the variables eliminated by 
the constraints are chosen so that the integrand is bounded in that region. 

Our second step then is to fix, by an integration in some c?e i? the closure 
constraint. We will need to be able to bound below the length of at least one 
edge of the chain. So define 

©g_ £ = { 9 € Q% } I abs(L(0) - b) > e } 

Q%;*> = {0ee£ } I abs(L(<?) - 6) < e, e, > (& - e)/n }, 

with i G {1,2, ...,n}, and e a small positive constant, < e < b, depending 
on b. Each chain in O^ has the property that its ith edge has length at 
least (b — e)/n. Any chain, with n edges and a total length differing from b by 
not more than e, must have such an edge. The sets Op* , i = l,2...n are 

not disjoint, but combine with ©p _ e to cover 0^ . Chains in Op _, will not 
contribute to the integral. It follows that 

Jn < E/ 5(b-L(§))5W(X k e k ) d ^---de n dxi 
~iJe£f eie 2 ...e„ 

< ^Ef, Ab-m) ^-^-^ ^. as) 

(0 - e) f— ' y e («.«) eie 2 . . . e_i . . . e„ 

where a — i subscript indicates that element is left out of a product or sum. 
Op is the set of polygons with a long ith edge (that is, the set of chains in 

0p'^ with x n+ i = x\). We have carried out the integral cfe^ 2 ) (Efce fe ) and 
used the bound on e,. 
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The third step is to eliminate an edge length parameter, using b — L(6), the 
length constraint. Let fa denote the angle made by edge to a fixed direction 
in the plane. In polar coordinates Equation ( |l3| ) is 

11 n f 

Jn < Ti nT/ 6(b - L(9)) deidfa . . .de-idfa-i . . .de n d(/) n dxi. (14) 

( b - £ ) Je^:> 



For the polygon to close 



eisin(fa) = -^e fc sin(</> fe ), (15) 

k=i 

n 

e,:cos(0i) = -^2e k cos(fa), (16) 



fe=i 

k^i 



and consequently 

n 

b - L(6) = b - ^2 e k - efc cos(^ fe 



fc=i 



Integrating dej for some j may lead to an unbounded integrand. In order to 



control this, we partition Q^'^ on its angle variables. Let 



{e i ^<\fa-fa\<^ } 

A polygon in Q^'^' has the property that the jth edge "turns back" from the 
direction of the long ith edge. There must be at least one such edge for the 
polygon to close. The sets Q^'^' , j = 1, 2 . . . n, j ^ i are not disjoint but their 
union covers Q$f. From Equation © 

n n n f 

Jn < 77 S{b — L{6)) deidfa. . .de-idfa-i .. .de n d(j) n dxx. 

We may now apply the integral dej to the delta-function S(b — L{6)). We 
transform from e, (f> to e', 0' where fa k = fa and e' fc = e k for 1 < k < n, k ^ j, 
and faj — <j)j and 

e'j — ej — ej cos(0j — 0i). 
The Jacobian of the full transformation e, <fi — > e', 0' is just 

J- 1 (e,0-,e',0') = 

= l-cos(fa - fa) -ejsm(fa - fa)-^-. (17) 

Repeated use of Equations ( |l5| ) and ([l6]) gives 

-sm(fa - fa) 
dej ei 
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in Equation ( |l7| ) and then using 7r/2 < \<pj — 4>i\ < 3ir/2, we have J^ 1 > 
1. The angle partition was needed to control this function. We can replace 
S(b — L(9)) dej by one, and restrict the integration domain to polygons of length 
b, ie set e = 0. We obtain the simplified bound 

J n < — / ^ de\d(j)i . . . de_jd<f)j . . . de-id<fi-i . . . de n d4> n dx\. (18) 

The last step is to bound the integral in Equation (p"8|). Enlarge 9^'^] to 
allow each variable to range independently over its full domain, keeping only the 
bound on total edge length, L(9) = b, and requiring x\ to remain in V. This will 
include polygons with crossing edges and allow the polygon to overlap the border 
of V. The integral dx\ gives a factor A(T>). Each angle variable ranges over 
to 2tt contributing (27r)™ _1 . The edge integrals are over the (n — 2)-dimensional 
tetrahedron 

e\ + e 2 + ■ ■ ■ e_j + . . . e_j + . . . + e n < b — b/n 

of volume less than b n ~ 2 J(n — 2)!. Combining these factors with a factor of 
(n — 1) from the sum over j, we obtain the bound on J n given in the Lemma. 
This is the end of the proof of the Lemma. 

Equation (0) is obtained by evaluating the integral over b in Equation ( [To| ) 
with the bound from Equation dl3), and the Theorem follows directly from 
Equation (||). 
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Figure 1: (A) A state x 01 the Arak process (B) The discontinuity set 7 of (A). 
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Figure 2: (A) A set of lines I intersecting V (B) an admissible graph drawn on 
the set I (C) one of the two colourings of V with discontinuity set given by the 
graph in (B). 




A B C 

Figure 3: Updates in the Markov Chain Monte Carlo. Dashed and solid edges 
are exchanged by the moves, which are reversible. (A) Interior vertex birth 
and death (B) move a vertex, and (C) rccolour a region by swapping a pair of 
edges. In an extra move, not shown, a small triangle may be created or deleted. 
Further move types are used to update boundary structures. 
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Figure 4: Binder parameter Ud (see text), regressed with cubic polynomials. 
Curves correspond to distinct box-side lengths d. The maximum likelihood fit, 
constrained to intersect at a point, is shown. Error bars in this and all other 
graphs are la. 
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Figure 5: The Binder parameter data of Figure |4| rescaled with Ising critical 
exponents. The regression is a cubic polynomial. X43-4 = 38.5 for the fit is 
acceptable. 
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Figure 6: The magnetisation mrf(T), regressed with cubic polynomials. 
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Figure 7: The magnetisation data of Figure |6| rescaled with Ising critical 
exponents. The regression is a quartic polynomial. The value of the \ 2 statistic 
shows that the fit is a poor one. 



20 



► 






~_ 4 

<■ . . . 







T=0.62 0.64 0.66 















W 





0.68 0.70 0.72 



Figure 8: A selection of states equilibrated in a box of side d = 12 at temper- 
atures below and above the estimated critical temperature T c ~ 0.6665(5). 
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